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Abstract 



The EPDifF equation governs geodesic flow on the diffcomorphisms with respect to a cho- 
sen metric, which is typically a Sobolev norm on the tangent space of vector fields. EPDiff 
admits a remarkable ansatz for its singular solutions, called "diffeons," whose momenta are 
supported on embedded subspaccs of the ambient space. Diffeons are true solitons for some 
choices of the norm. The diffeon solution ansatz is a momentum map. Consequently, the 
diffeons evolve according to canonical Hamiltonian equations. We examine diffeon solu- 
tions on Einstein spaces that are "mostly" symmetric, i.e., whose quotient by a subgroup of 
the isometry group is 1-dimensional. An example is the two-sphere, whose isometry group 
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S0(3) contains S^. In this situation, the singular diffeons (called "Puckons") are supported 
on latitudes ("girdles") of the sphere. For this symmetry of the two-sphere, the canonical 
Hamiltonian dynamics for Puckons reduces from integral partial differential equations to a 
dynamical system of ordinary differential equations for their colatitudes. Explicit examples 
are computed numerically for the motion and interaction of the Puckons on the sphere with 
respect to the norm. We analyse this case and several other 2-dimensional examples. 
From consideration of these 2-dimensional spaces, we outline the theory for reduction of 
diffeons on a general manifold possessing a metric equivalent to the warped product of the 
line with the bi-invariant metric of a Lie group. 
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1 Introduction 

1.1 Motivation and Problem Statement 

When the Schumacher-Levi comet broke into a string of fragments that collided with Jupiter 
in July 1994, each of the impacts created a circularly symmetric expanding "ripple" spreading 
out into the gaseous sphere of Jupiter. If the collision had occurred instead with an Earth- 
like planet which was entirely covered by a shallow layer of water, it would have produced 
a sequence of strongly nonlinear circular traveling-wave ripples spreading out from a point 
source on the spherical surface, propagating to the antipodal point, then reflecting back in the 
opposite direction and colliding head-on with the trailing waves. 

Neglecting planetary rotation and assuming the energy was mainly kinetic in this situation 
would produce an approximate description of these strongly nonlinear shallow water waves as 
solitary waves governed by the Camassa-Holm (CH) equation [3], written on the sphere. The 
momenta for these solitary waves on the surface of the sphere would be concentrated on a series 
of coaxial circles which could each be regarded as a latitude encircling the axis running from 
the point of impact to the antipodal point. These concentrations of momentum on latitudes are 
singular solutions of the EPDiff equation, which is a geometrical extension of the CH equation 
for shallow water waves to higher dimensions. 

Previously, the CH equation and its extension EPDiff have been solved analytically for the 
interactions of such waves on the real line in one dimension [3], as radially-symmetric rotating 
concentric circles on the two-dimensional plane [8] and numerically in two dimensions and three 
dimensions, for a variety of initial value problems [10]. In cases with one-dimensional linear, 
or radial symmetry, these dynamics reduce to ordinary differential equations (ODEs). Here we 
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consider the corresponding reduction for singular wave motion on the surface of the sphere for 
EPDiff with respect to the norm and we discuss its generalisations for other surfaces of 
constant curvature. 

In each case, we derive the symmetry which reduces the dynamics of the singular solutions 
of the EPDiff equation to a system of ODEs. We also provide explicit numerical results for 
the interactions of these singular EPDiff solutions in the case of Puckon motion as concentric 
latitudes moving on the sphere. 



1.2 Review and Generalisation of the Camassa-Holm Equations to EPDiff 

The dispersionless Camassa-Holm (CH) equations [3] for one-dimensional shallow water waves 
arise as stationary points of the kinetic energy functional given by the Sobolev (2,l)-norm 

i[u] = / {uix)"^ + Ur.ix)'^) dx , 
Jr 

subject to velocity variations in the Euler-Poincare form [6] 

6u = C+ [C, u] , 

where ^ S g, and g consists of the vector fields on the real line, with Lie bracket [•, •] given 

by 

[v, w] = vwx — wvx = ad^w . 

The dispersionless CH equations themselves are thus given by the following system for velocity 
u and its dual momentum m, 

m = — ad*m = —{mdx + dxm)u , 

where m = — = u — u^x = u + An , 
ou 

u = G*m= I G{x — y)m{y)dy , 
Jr 

and G{x) = ^e"'^' is the Green's function for the Helmholtz operator 1 — on the real line. 
We remark that here we are using the convention that the Laplacian in ID is given by 

Au = —Uxx- 

The sign comes from regarding the Laplacian as d*dx, where d* is the adjoint of the 
derivative operator dx- 

To put the dispersionless CH equations onto a general Riemannian manifold (M, (•,•)) with 
dimM = n and Levi-Civita connection V, we work with the functional on the space of weakly 
differentiable square-integrable vector fields 

l[u] = [ + |Vup)dfo/. 

Jm 
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Again, we consider the stationary points of H. with respect to variations of the Euler-Poincare 
form 

5u = C+ [C, u] . 

This requirement yields the EPDiff" equations, for "Euler-Poincare equations on the diffeo- 
morphisms," given by the system, 

m = — ad*m , where m = — = u + A^u . (1) 

5u 

Here we denote 

a'^ = v*v 

for the connection Laplacian with respect to the metric and we assume homogeneous boundary 
conditions. 



Diffeons: Singular solutions for EPDiff. The EPDiff equation (1) has Lie-Poisson Hamil- 
tonian structure in the momentum variable m and satisfies a remarkable ansatz for its singular 
solutions, which are called "diffeons," 

N 

m{t,x) = \2 P^{t,s^)6{x - Riit,Si))dsi. (2) 

The diffeon singular solutions of EPDiff are vector-valued functions supported in M on a set 
of A'^ surfaces (or curves) of codimension (n — k) for s G M'^ C M with dimM'^ = k < n. 
For example, diffeons may be supported on sets of points {k = 0), curves {k = 1), or surfaces 
(k = 2) in three dimensions. These support sets of m move with the fluid velocity u = G * m; 
so the coordinates s £ C M are Lagrangian fluid labels. The Green's function G for 
the operator 1 -|- is continuous, but it has a jump in its derivative on the support set that 
advects with the velocity u = G*m under the evolution of EPDiff. In fluid dynamics, a jump in 
derivative which moves with the flow is called a "contact discontinuity" [12]. The relationship 
of the singular solutions of EPDiff equation (1) to fluid dynamics is illuminated by rewriting 
equation (1) in "Riemann-invariant form," 

— (m ■ dx (8> dvol) = , along — = u = G * m . 
dV ' dt 



The diffeon singular solution ansatz (2) was first discovered as the "peakon" solutions for CH 
motion on the real line in [3]. This was generalised to motion in higher dimensions in [9] and 
was shown to be a momentum map in [5]. As a result of the singular solution ansatz being a 
momentum map, the 2A^ variables Pi,Ri, satisfy Hamilton's canonical equations. In general, 
these are integro-partial differential equations for the canonically conjugate diffeon parameters 
Pi{t,Si),Riit,Si) in (2) with i = l,...,N. 
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1.3 Aim of the paper 



This paper will be concerned with examining diffeon solutions of the EPDiff equations (1) on 
Einstein manifolds that are "mostly" symmetric, i.e., that have a group acting by isometries 
so that the orbits have co-dimension 1 in the manifold except on a set of measure zero, and 
hence the quotient of the manifold by the group is 1-dimensional. 

Thus, we are interested in identifying and analysing cases where imposing an additional trans- 
lation symmetry on the solution reduces the canonical Hamiltonian dynamics of the singular 
solutions of the EPDiff on Einstein manifolds from integro-partial differential equations to or- 
dinary differential equations for Pi{t), Ri{t). We shall then analyse some of the properties of 
those canonical equations and examine their numerical solutions. 

We begin by noting some simplifications of the EPDiff equations when restricted to Einstein 
spaces. These simplifications take advantage of the relation between the connection Laplacian 
and the Hodge Laplacian, the latter of which is easier to compute. Knowing the EPDiff 
equations in 1-dimension, we move up to the Einstein spaces in two dimensions. 

The first manifold we study is the 2-sphere, which has the familiar group action formed by 
rotations about a fixed axis. Thus we construct "Puckons," which are singular solutions of 
the EPDiff equations supported on concentric circular latitudes of the 2-sphere. The name 
"Puckon" (as opposed to "peakon" ) arises from the famous boast by Shakespeare's character, 
Puck, in A Midsummer Night's Dream, that he would, "put a girdle round about the earth 
in forty minutes."^ The case of rotationally symmetric singular solutions of EPDiff in the 
Euclidean plane has already been discussed in [8], so we proceed to examine the hyperbolic 
plane, which has a rich isometry group. However, the diffeon dynamics on the unbounded 
hyperbolic plane affords little opportunity for multiple diffeon interactions and does not admit 
periodic behaviour. Richer opportunities for diffeon dynamics and interactions in hyperbolic 
spaces are offered in the context of Teichmiiller theory, which is, however, beyond the scope of 
the present work. 

From consideration of these 2-dimensional spaces, we sketch how one might develop the theory 
for the general manifold which possesses a metric equivalent to the warped product of the line 
with the bi-invariant metric of a Lie group. 

1.4 Plan of the paper 

After briefly recalling the essentials needed from the theory of Einstein manifolds in section 2, 
we begin in section 3 by reducing the singular solutions of EPDiff to a canonical Hamiltonian 
dynamical system for the simplest Einstein manifold - namely, motion on a sphere of singular 
EPDiff solutions supported on concentric circular latitudes, or "girdles." These new singular 
solutions girdling the sphere are the "Puckons." Although the canonical reduction is guaranteed 
by the momentum map property of EPDiff, the reduction of its singular diffeon solutions for 

^Thanks to J. D. Gibbon for reminding us of this quote and suggesting the name, "Puckon" for these solutions. 
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Puckons is given in detail, so it may be used to confirm the numerical results for the interactions 
of Puckons on the sphere described in Section 2. Section 5 generalises the Puckons to other 
surfaces that are Einstein manifolds with a translation symmetry. Section 6 incorporates these 
ideas into the theory of warped products. 



2 The EPDifF Equations on Einstein Spaces 

The operator 

+1 

arises in the following context in Riemannian geometry. 

Definition 2.1 On any given Riemannian manifold, the musical isomorphisms are defined to 
be the Riesz representation and inverse maps with respect to the metric: 

(j : T*M — > TM 

b : TM — > T*M 
v\,{w) = {v,w). 



In Riemannian geometry, the Levi-Civita connection respects the musical isomorphisms, namely 
for v,w £ X{M) and a £ n^{M): 

(V,a)« = V,(a«). 

Thus, the musical isomorphisms identify vector fields with 1-forms and allow them to be dif- 
ferentiated in the same way 

The connection Laplacian on 1-forms satisfies the Bochner-Weitzenbock formula. 



Theorem 2.2 (Bochner-Weitzenbock) On a Riemannian manifold with Levi-Civita con- 
nection V, any 1-form a satisfies 

A'^a = A'^a + Ric(a) 

where A'^ = d*d -\- dd* is the Hodge Laplacian formed from the exterior derivative d on forms 
and Ric is the Ricci curvature operator. 



Now we restrict our attention to special forms of manifold, namely Einstein Manifolds. 
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Definition 2.3 An Einstein manifold is a Riemannian Manifold which satisfies 

Ric = kt 

for some constant k. 



On an Einstein manifold, one may scale the metric so that k can be replaced by —1,0,1 
depending on whether k is negative, zero or positive respectively. Thus on an Einstein manifold, 
the Bochner-Weitzenbock formula becomes 

= + sign(A:)l. 

We restrict our study further to Einstein manifolds with positive k (we shall call these positive 
Einstein manifolds) and scale to A; = 1. Thus we have 

A'^ = A^ + 1. 

The implication of this for the EPDiff Lagrangian on vector fields is as follows, for homogeneous 
boundary conditions: 

e[u] = [ + iV-up) dvoZ 

= / (|ub|2 + |VH|^)dt;oZ 
Jm 

(lit,, U[, + A^uj,) dvol 
(^ui,, A^ii[,^ dvol 

( = [ {\dUi,\'^ + \d*Ui,\'^)dvol). 

Hence, finding the stationary points of Lagrangian I with the usual Euler-Poincare constraints 
on the variations implies the EPDiff equation system of equations (1) in the form 

m = — ad^m , where m = l^^u\, . (3) 

These equations generalize the EPDiff equations to Einstein manifolds. 





3 The EPDiff equations on the Sphere 

In two dimensions, the only manifold with positive Einstein constant is the standard round 
sphere which we shall regard as the Riemann sphere. We use stereographic projections to 
identify the complex plane with the sphere whose North Pole is removed. This is equivalent to 
putting the metric 

4 / 2 2\ 4 

3 = TT^ — 2 I 2\2 (^^ + ) = ——-^dzdz 

(i + x^ + y^)^ ' \\^\zYY 

on the plane. 
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3.1 Rotationally Invariant Solutions 



We shall examine vector-field solutions u of EPDiff in (3) witfi 



AV,/3)= / (3 



S2 ^ ' JC 

for any smooth 1-form (3 where C is a circle of latitude on the sphere. The quantity m = A^U|, 
is a distribution, defined by its integration against a smooth function. Thus, we seek weak, or 
singular, solutions of EPDiff on the sphere. We change coordinates on the sphere to assist us 
in this search. The metric on ]R^\0 in polar coordinates is dr^ + r'^dO. Thus we can regard the 
sphere minus both poles as ]R^\0 with the metric 

gs = (dr^ + r^d6'^) , where p = ^ . 

\ / M _|_ 



Green's function for the Helmholtz operator on the Riemann sphere. We now seek 
solutions to the equation 

A^itb = S{r - R) {kidr + fead^) (4) 

where R, ki, k2 are constants. Let us assume that the velocity one-form u\) is invariant under 
rotations of the sphere about the axis joining North and South Poles and is radial, i.e. we must 
solve 



U\, = a(r)dr, (5) 
Let us first solve the Green's function equation 



A<i(G(r,i?)dr) = ^i^— ^dr, 
r 

for one recognises that the single diffeon velocity is proportional to the Green's function, i.e., 

a(r) = PG{r,R). 

To solve explicitly for the Green's function in our present case, we begin by recalling 

or \ Ar or J 

Integrating the Green's function equation using this expression yields 

G(r, R)dr = (— + + — — ^ , ) dr , 

\ r r(l + r^j r(l + max(r, /tj^ j y 

for constants Ai,A2. In particular, we can remove the singularities at r = and r = oo by 
setting Ai = and 

A.- ' 



R{1 + i?2) ' 
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whence 

Similarly, the solution to 

A<^G(r, R)de = R6{r - R)de 
for the angular diffeon velocity component is 

Both G and G are continuous over the sphere, but each has a jump in derivative at r = i?. 
Thus, we have proved the following. 

Proposition 3.1 (Radial Green's function on the sphere) The solution to (5,6) for the 
radial Green's function on the sphere is 

2Pmin(r,ii)2 

Solution ansatz for EPDiff on the sphere. Following [9], we propose a solution ansatz 
for EPDiff velocity on the sphere as the following superposition of Green's functions, 

N 

r? = ^P,G«,dr (7) 

i=l 

where we denote 

GR{r) = RG{r,R). 

and Pi, Ri are 2N functions of time. The corresponding vector field dual to rj is 

^ 1 

u = y^^PiGR^dr (8) 

We pair the arbitrary smooth vector field w = fdj- + gdg on with the EPDiff equation using 
the solution ansatz given by the one- form rj in equation (7). A direct calculation yields the 
following equations for the diffeon parameters. 

Proposition 3.2 (Diffeon parameter equations on the sphere) 

The equations for diffeon parameter evolution on the sphere are: 

^ f P- \ 

= T.[piR,)^piR^)2R^^^''^'^^^) (9) 

+ 2R,p{R,)G{R^,Rj)\ . (10) 



^^piR^ypiR^y (sr^^"'^^) 



r=R, 
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Proof (by direct calculation): 

= ^^(^j^ A%{w)dvol^ + A'^ri{[w,u])dvol 



d_ 
dt 



E 



JO 



5{r - Ri)f{r,e)Pip^rdrde 



1-K ^ t-oa 



JO 



N 



f{r, 0)dr + g{r, d)de, J2 P'^PjGr^ (O^r- 



o^rdrd9 



+ Y.^' [PjG'ji^{Ri)fiRi,0) + 2PjRip{R,)GR^{Ri)f{R„d) - P,GR^{R,)fr{Ri,e)) P.dO 

i,j=l 

= V/ {fr{R^,e)R^P^ + P^f{R^,9))d9 

N ^2n 

+ Y [^^^'r. iRi)f{Ri^^) + 2PjRip{R,)GR^ {R^)f{R^, e) - PjGr^ {R^)fr{Ri, 9)) P.dO 



Comparing coefficients of / and implies 



N 



RiPi — ^ 



^ \p{R,fp{Rj) 



jG{Ri, Rj 



N 



> piR,yp{R,)m, 



G'R(R^) + 2RipiRi)GR^{R^] 



This finishes the calculation of the diffeon parameter evolution equations (9,10). 







Proposition 3.3 (Canonical Hamiltonian form of diffeon parameter equations) 

Evolution equations (9,10) for the diffeon parameters are equivalent to Hamilton's canonical 
equations, 



vr dRi ' vr dPi 



:ii) 



with Hamiltonian function ofP = {Pi, . . . , Pn) and R = (i?i, . . . , Pjy) given by. 



2 P{R^YP{R,? ■ 



(12) 
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As explained in [5], this reduction of EPDiff to canonical Hamiltonian form for the diffeons 
is guaranteed, since the singular solution ansatz (2) is a momentum map. However, we shall 
require the explicit results for the numerical solutions in section 4- 



Proof: We Legendre transform the Lagrangian i into the Hamiltonian H by setting 

H{m) = / Tn{u)dvol — £[u] 
Js2 



whence we find that 



with 



H{m) = - f m{u)dvol 
2 Js2 



m = A'^uu. 



We then express the velocity as the superposition of Green's functions, 

1 2P,mm{r,R,f 
^^p{R^rrR,{l + r^){l + R^^■ 

We set ki = Pi / {p{Rif Ri) , for i = 1, . . . , A^, and we evaluate the Hamiltonian on this solution 
with m = A'^r] as 

H{m) = - [ {r],m)dvol 

= - / / > ( — \, dr,A'^—^ — ^ \ dr ) p^rdr A de 
2 70 Jo ^j^,\ril + r^)il + Rh r(l + r2)(l + i?2) 

= o / / E C'^^R^ (^)d^' - ^M'') /''^dr A de 

2 Jo Jo ,^^.^1 

= o / ^hGn,{R,)k,R,de 
^ vr ^ P,PjG{Ri,Rj) 

2 p{R.rp{Rjy • 

The canonical equations for this Hamiltonian now recover the diffeon parameter evolution 
equations (9,10). | ^/ 1 



Remark 3.4 (Remark on constant factors) For convenience and brevity of writing, in 
what follows we shall often absorb into the time units any constant factors - such as vr in 
equations (11-12) - arising from integration over symmetry directions. 
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Definition 3.5 (A^-Puckon) The singular solution of EPDiff 

^^AS = -ad;,A^, (13) 
is given on the Riemann sphere by a vector field r] satisfying 

^ P 

= E^(#^^(^-^Od.. (14) 

The support set of this (weak) solution of EPDiff on the Riemann sphere is a set of circular 
latitudes (girdles) at radii Ri{t) with conjugate radial momenta Pi{t), where i = 1,...,N. 
Equation (I4) constitutes a solution ansatz for the velocity vector field ij, which will be called 
an N -Puckon solution. 



3.2 The Basic Irrotational Puckon 



Purpose. This subsection derives explicit solutions for the parameters P and R of the single 
Puckon without rotation and examines its behaviour upon collapsing to the poles of the sphere. 



The single irrotational Puckon. Let us consider the case where = 1 and examine the 
motion of the basic Puckon without rotation. For = 1, the Hamiltonian (12) is given by 

1 P'^G{R,R) _ „2n2 

which follows because the Green's function in this case is 

mm{R,R)'^ 2 



G{R,R) = 2- 



'R^{l + R^)^ (l + i?2)2- 
Thus, upon restricting ourselves to the constant level set of the Hamiltonian defined by 

H{P, R) = — 

for constant K, we may solve for the momentum variable, 

2K 



P = Kp{R) 



(1 + R'') ■ 

Next, we note that the canonical coordinate equation 
integrates to 

R = tan \ —t + -] . (15) 
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Consequently, the single Puckon momentum is found as 

P = K + Kcos{Kt + e). (16) 

So far we have been using the stereographic projection to provide charts for the sphere. How- 
ever we may easily pass from the coordinates (r, 6) obtained from stereographic projection to 
latitudinal-longitudinal coordinates (0, 6) where r and (p are related by 

r = tan 

By this token, the colatitude (j) = ^{t) of the peak of the Puckon evolves linearly in time as 

^{t) = Kt + £. 




Proposition 3.6 (Irrotational Puckon Solution) The velocity vector field v = rf' gener- 
ated by the irrotational Puckon motion is given in stereographic coordinates by 

P min(r, R) 



u 



-{l+r')dr. 



4 max(r, R) 

This appears in latitudinal-longitudinal coordinates on the Riemann sphere, as follows, 



u 



K{l + cos{Kt + e)) ™™ 


( 


tan 


(t) 




tan (ft +1)1 


) 6 


^ max ^ 


tan 


(1) 




tan (ft + 1)1 





Proof: This result is obtained by direct substitution of solutions (15) and (16) for the diffeon 
parameters into the solution ansatz (8). 

Remark 3.7 (Puckon peak) The peak of the Puckon occurs when 

(f) = $(t) = Kt + e, 
(modulo behaviour at the poles), whence 

, , \K\(l + cos(Kt + e)) 

\u\a2 = . 

I lb 2 

We still need to examine precisely what happens to the dynamics of a Puckon as it collides 
with itself at the poles. 

Proposition 3.8 (Puckon behavior at the poles) The Puckon bounces elastically, when 
it collides with itself at the poles. 

Proof: We notice that R = K/2 / at i? = and by setting U = l/R we see that at C/ = 0, 
U = K/2 7^ 0. Similarly, one can show that P = at both poles. Thus, the Puckon velocity 
is finite at the poles. Since Hamiltonian motion is time-reversible, the Puckon must bounce 
elastically, when it collides with itself at the poles. 
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3.3 Rotating Puckons 



So far we have concentrated on singular EPDiff solutions on the sphere moving with only a 
radial component. Now we turn to examine the rotating A^-Puckon, i.e. the solution of the 
following system of equations for tj, cf. equations (4), 

= -ad;,A''>,, (18) 
in which Mi(t) with i = 1, . . . , N , are the angular momenta of the Puckons. 

Proposition 3.9 (Canonical Hamiltonian equations for the rotating A^-Puckon) 

The equations for the rotating N-Puckon may be expressed as: 

Ri = Y.PoFARi)R^R^ = ^,^ (19) 

^ dH 
Pi = -Yl {PiPiF',{Ri)R3R^ + PiPjFi iRi)Rj + MiM,F^{R,)) = - ^ , (20) 

Mi = 0, (21) 
where Fi{Rj) = Fj{Ri) is defined in terms of the Greens function G as 



Fi{Rj 



G{Ri,Rj) 



p{RifRiP^{Rj)Rj ■ 

The last equation (21) expresses the conservation of the angular momentum of each Puckon. 
With Mi constant, the first two equations (19) and (20) comprise Hamilton's canonical equa- 
tions with symmetry-reduced Hamiltonian 



1 ^ 

F(P,R;M) = -J2iP^PJR^Rj + M^Mj)Fi{Rj). (22) 
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Proof. First, we note that the velocity one- form corresponding to the momentum singular 
solution ansatz (17) is given by, 

" G(r,R,) 



Thus the velocity vector field u dual to rj is given by 

N 
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To make our calculations more transparent, we collect terms and define notation as 

G{r,Ri) 



so that 



and 



N 

u = Y,Fi{r) {R^P^r^r + Mi^e) 



i=l 

N 



77 = ^ Fi{r)p^r {RiPidr + MirdO) . 

i=l 

Note that the expression for Fi is symmetric in r and Ri, which implies the permutation 
symmetry, 

Fi{Rj) = Fj{Ri). 

Now let w = fdr + gdg be any smooth vector field on S^. Then, one computes the pairing 

/ ^''r]{w)dvol = V/ {Pif{Rue)+Mig{R,,e))de. 
Hence, the left hand side of the EPDiff equation becomes 

^ (^^^ AMw)dvol^ = J2 fj {PifiRi, 0) + P^Rifr{Ri,6) + IvUg{Ri, 9) + M,Rigr{Ri, 6)) d9. 

Now we need to calculate its right hand side. 

A'^r]{[w, u])dvol. 



For this, we write the commutator, 

N 

[w, u] = {F^R^P^rf + F^RiP^f " F.RiPirfr - FiMJe) 5. 

i=l 
N 

+ if^^^i - P^PiPi^dr - FiM.ge) de. 
1=1 

Now, by Stokes' theorem, when we integrate over the whole sphere, the contribution due to 
f9,ge will be zero because f,g, the only components of the integration to depend on 9, satisfy 
/(r,0) = /(r,27r) etc. Thus 

d 

2 



A ri{[w, u])dvol 

V r {P^F'{R,)RjP,Rif{R,, 9) + P,Fj{R,)RjPjf{Ri, 9) - PiFj{R,)RjPjR,fr{R,, 9)) d9 
+ V / {M,F'{Ri)Mjf{Ri, 9) - MiFj{R,)RjPjR,griRi, 9)) d9. 
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Now to find the evolution equations for Pi,Ri,Mi we need only compare the coefficients of 
f,fr,g,9r occurring in 

= ^(^^ A'^7]{w)dvo?j + A%{[w,u])dvol. 

From this comparison, one finds 

N 

= Pi + Yl {PiF■{Ri)RJPjR^ + P^Fj{R,)RJPj + MiF'^{Ri)M^ , 

n 

= PiRi — ^ (^PiFj{Ri)RjPjRi^ , 
= Mi, 

N 

= MiR,-Y,{M,Fj{Ri)RjPjR^, 
i=i 

in which the second and fourth equations provide the same information. We therefore obtain 
the evolution equations (19-21) for Pi,Ri,Mi given in the statement of Proposition 3.9. | ^y | 

Remark 3.10 (Verifying the Hamiltonian) A simple check shows that if the Hamiltonian 

H{m) = / m{u)dvol — i[u] 

is the Legendre transform of the Lagrangian I then 

H{m) = (^A'^ri,ri^ dvol 

= 27rF(P,R;M) , 

which is the Hamiltonian for the rotating N -Puckon in formula (22) of Proposition 3.9. As 
explained in [5], the reduction to canonical Hamiltonian form is guaranteed, since the singular 
solution ansatz (2) is a momentum map. 



3.4 The Basic Rotating Puckon 

Proposition 3.11 (Extremal radii of the basic rotating Puckon) 

The motion of the basic rotating Puckon lies between the maximum and minimum values of R 
on the Riemann plane given by, 

2K ± V4:K'^ - M2 

-'^max/min i {■^'^J 

assuming that 2K > M > 0. 



17 



Proof: We = 1 in Hamiltonian (22) to find 

H{P, R) = -^{1 + R^)\P^R^ + m2). 

We notice that for a rotating Puckon with H = = const, the girdle of the Puckon cannot 
have zero radius unless M = 0, in which case we return to the irrotational Puckon. Likewise, 
the girdle radius cannot be infinite unless again M = 0. Thus, a rotating Puckon with M 7^ 
is constrained to lie between a maximum and a minimum radius. To find these extremal radii, 
we must solve 

dH P, 

H = -^{1 + R^fiP'^R^ + M'^) = K^. 
The first equation can only be solved by P = 0; for which the second becomes 



[\^Ry = K\ 



16R'^ 

Assuming that K and M are both positive, we find that the maximum and minimum values 
of R are the roots, 



_ 2K± V4:K'^ - M2 

^max I min 

This proves the proposition. | ^/ 1 

Remark 3.12 (A single critical point) We observe that if M = 2K, then the Puckon is 
radially static at the equator and P is identically zero. The solution {P,R) = (0, 1) is the only 
critical point of the Hamiltonian H unless M = in which case the set defined by P = is the 
critical manifold. 

Proposition 3.13 (Periodic motion of the basic rotating Puckon) 

The motion of the rotating Puckon is periodic, with period I'k I\J2H determined by the constant 
value of the Hamiltonian, H . 



Proof using the colatitude representation of periodic rotating Puckon motion: 

To investigate the motion of the rotating Puckon, we pass from stereographic coordinates 
to longitude-colatitude coordinates in which we put 

R = tan — . 
2 

As we will see later in the general case, this does not alter the situation a great deal. For 
example, the momentum is given later in equation (32), after substituting sin for ^. The 
Hamiltonian for the rotating Puckon on the sphere in longitude-colatitude coordinates is 

H{P,^) = Up^ + ^), (24) 
2 V sm 5> / 
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in which <I> and P are canonicahy conjugate variables. Along the motion of the Puckon, H is 
constant, say H = K^. Thus, taking the positive branch for P 



sin^ $ 

yields the equation for the colatitude 



M2 

$ = P = W2K2 



sin"' $ 

Integration of this ODE implies that 

sin$d$ 

t 



<i>{0) V2A'2 sin2 $ - M2 

1 / ( K^/2cos^{t)\ /a'V2cos^>(0)\\ 

= ■= arccos — =■ — arccos — =- 

Ky/2 \ \ V2K2 _ J \ V2K2 - A'P J J 

Thus the colatitude for the single rotating Puckon is given as a function of time by 



/2i^2 cos2$(0) /2A'2sin2$(0) -M2 

^(0 = V 2i^2_M2 -osiKtV-2) - y -(E:tV2). 

These rotating Puckon solutions are periodic with period 2tt/{K\/^), determined by the con- 
stant value of the Hamiltonian. | ^/ 1 

3.5 Puckons and Geodesies 

Proposition 3.14 (Geodesic motion of a point on the girdle of the rotating Puckon) 

A point on the girdle of the rotating Puckon moves along a great circle at constant speed equal 
to K\f2 determined by the value H = of the Hamiltonian, H. The normal to its plane of 
motion is inclined to the South Pole at the angle vr — arctan(2i2max) with Rmax given in (23). 



Proof: For a rotationally symmetric surface M with metric 

g = dr'^ + ^{rfde, 

the equations for a curve x = {R{t), @{t)) to be a geodesic are equivalent to 

^2+^202 ^ ^ (25) 

^29 = const. (26) 
So in the case of the sphere with metric g = dcj)"^ + sin^ (f>d6, if we define 

M2 

e = uei'^) = -^ (27) 
sm 9 



19 



then we obtain a curve on the sphere with coordinates x = ($, 0) and tangent vector 



x{t) = (6,9). 



Now, the speed is given by 

|xp = g{x,x) 



6^ + (sin^ 

p2 



sin^ $ 

= 2H{P, ^>) = constant = 2K'^ . 

Thus a point on the geodesic has constant speed equal to \/2H and by (27) we have 

(sin^ <I>)6 = = constant. 

Consequently, x/{K\/2) satisfies equations (25) and (26), thereby determining a geodesic on 
the sphere. This means that a point on the girdle of the Puckon moves along a great circle at 
constant speed equal to \/2H and the normal to its plane of motion is inclined to the South 
Pole at the angle vr — arctan(2i2maa;)) with Rmax given in equation (23). | ^/ 1 



3.6 Further Hamiltonian Aspects of Radial Solutions of EPDiff on the Rie- 
mann Sphere 

Proposition 3.15 (Lie-Poisson Hamiltonian form of EPDiff on the Riemann sphere) 

The radially symmetric solutions of EPDiff on the Riemann sphere may be written as 



— ( \ = v( ^' \=v( ^^/^"^ 
dt\rme J \^ J V SH/6{rme 



^( )=^( ) > m 



where T> is the skew- symmetric Hamiltonian operator given by 

X> = f ~ {"^rdr + drnir - 2mrpr + ^) -rmgdr \ ^29) 
y (y—drrruQ + 2mQr'^p — mo) J ' 

and the velocities Ur and uq jr are given by the variational derivatives of the Hamiltonian, 

5H ue 5H 



Ur 



5mr ' r 6rmg 



Equations (28,29) provide the Lie-Poisson Hamiltonian form of the EPDiff equation for radially 
symmetric dynamics on the Riemann sphere. 



20 



Proof: In the case that 



u 



G{r,Ri) 
^ p{RiYRip^r 

UrOr H 00. 



{RiP^rdr + Mide) 



where the Pi and Ri satisfy (19) and (20), the associated momentum density is 
where 



N 



rUr 



N 



Sjr - R{) 
^ p{R^?Ri 



Let wi = -^dr and W2 = -^d^ be two vector fields on S^. These satisfy the commutator 
relations, 



[wi,u] 



-ir-Or,UrOr H Og 

fu'r _ Urdrf _ fUr fu^ _ Updef 
1 1 22 22 



\dr+ (^-^]de, 



and 



[W2,U\ 



-^de,UrOr H O0 

UrdrQ gUr _ QUr uedeg\ 



p2y 



p p^r"^ p^r"^ J 



Thus, the EPDiff equations become, for wi, 
= f ( -^rniwi) + rn(\wi,u\)\ dvol 



2n noo 



JO 
2-K roo 



dnir f ^ ^ //m^ Urdrf ^fUr ^ fUr Updef 

dt p^r ^ \ p^r p^r p p^r"^ p^r"^ 



rmg 



fu'g _ fug \ 



p rdrd^ 



JQ 

—Qff + mr{ fu[. - Urdrf - 2furpr + — ^ ) ) drd^ 
JO \ r-^ 



From this we obtain the radial equation 
drrir 



dt 



(mrdr + drTHr — 2mrpr -\ -] Ur — rmgdr — . 

\ r J r 



p rdrd0 



(30) 
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Similarly, for 1^2, 

= j (^■^m{w2) + m{[w2,u])^ dvol 

Jo Jo \ P^'' V P p"^^"^ P^^^ / / 
= j j [ Qi 9 - rme [urOrg + 2gurrp — j I drdS 

Hence we have the azimuthal equation 

^ = i-d.rme + 2myp-me)u.. (31) 
ot ' 

Equations (30) and (31) provide the Lie-Poisson form of EPDiff on the Riemann Sphere. | | 

4 Numerical Solutions for EPDifF on the Sphere 
4.1 Overview 

We present numerical solutions to both the EPDiff partial differential equations (28,29) and 
the corresponding ordinary differential equations (19,20) for Puckons. Instead of using the 
stereographic projection, the numerical solutions were calculated on the sphere with colatitude- 
longitude (<?!>, ^) coordinates (and canonical variables instead of R,P). The equations in 
{(j), 9) coordinates on the sphere are obtained from those in Section 5 using ip = sin, so that 
g = d(/>2 + sin2(/)d6'2. 

We also introduce a length scale a for our numerical solutions by effectively changing the radius 
of the sphere from 1 to a. In this case a^A'^ = 1 + a^A^ so that the Green's function is 



G'(<^,^) = f 



tan I cot f ) , (/)< ^> 
f cot| 



tan tcoti]^^°' , (b> <^ 



4.2 Numerical specifications. 

Numerical simulations for diffeons on the sphere were performed using both PDEs and ODEs. 
In simulating the Lie-Poisson partial differential equations (28,29) in colatitude-longitude coor- 
dinates in Proposition 3.15, fourth-order finite differences were used to calculate spatial deriva- 
tives, and the momentum was advanced in time using a fourth-order Runge-Kutta scheme. The 
Hamiltonian was verified to be conserved to within 0.01% of its initial value for the simulations 
with smooth initial velocity distributions and to within 5-10% for the simulations when using 
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D pi/3 2*pi/3 pi 

Figure 1: An initial Gaussian breaks up into Puckons, two of which are evident for the times 
shown. 

one or two Puckons as the initial velocity distribution. The ordinary differential equations 
(19,20) for the canonical variables in colatitude-longitude coordinates were advanced in time 
using a fourth-order Runge-Kutta scheme. The Hamiltonian was conserved for these ODE 
simulations to within 10~^%. For PDE simulations, the results will be shown as a velocity 
distribution changing in time, whereas the results of the ODE simulations will be shown as 
plots of the time evolution of the canonical variables $,P. The length scale a was set to 0.1 
unless otherwise noted. 

Irrotational Puckons. We first consider irrotational Puckons {uq = 0). Figure 1 shows the 
evolution when the initial meridional velocity u^(0^ 0) is a Gaussian. As time elapses (t = is 
at the bottom of the figure, and later times are shown above it), a Puckon emerges from the 
Gaussian in Figure 1, and a second Puckon begins to emerge. In agreement with the prediction 
of equation (24), the first Puckon retains its height (and thus retains its velocity and canonical 
momentum P) as it approaches the South Pole. 

A rear-end collision between two irrotational Puckons is shown in Figure 2. As the plot of 
colatitudes shows, these two Puckons do not pass through each other. Instead, they bounce 
and exchange momenta. However, their momenta are not exchanged exactly, as the plot of 
P(t) shows. This contrasts with 2-soliton collisions for completely integrable PDEs in which 
momentum is exactly exchanged. 

A head-on collision between two irrotational Puckons is shown in Figure 3. In the PDE simu- 
lation, a vertical slope appears to form in finite time. Note that the Puckon velocities (i.e., the 
Puckon heights in the PDE simulation or the slopes of ^>(t) in the ODE simulation) remain 
finite and actually decrease to zero, whereas the equal and opposite canonical radial Puckon 
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time time 



Figure 2: Two Puckons undergo a rear-end collision. The PDE simulation agrees with the 
colatitudes ^ of the ODE simulation to within 1% and with the canonical momenta P to 
within 3%. 
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4 8 12 4 8 12 

time time 



Figure 3: Two Puckons undergo a head-on collision. 
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momenta diverge as the collision takes place. 



Rotating Puckons. Now we consider rotating Puckons (ug 7^ 0). Figure 4 shows the PDE 
evolution when the initial meridional velocity is zero and the initial azimuthal velocity uq 
is a Gaussian. Two rotating Puckons have emerged from the Gaussian by the time shown, and 
they are moving toward opposite poles. The Puckons each have small but nonzero azimuthal 
velocities. 

Finally, Figure 5 shows an ODE simulation of a single rotating Puckon when the length scale 
is a = 1 and the canonical variables are initially $ = 1.5, P = — 1, M = — 2. As discussed in 
Subsection 3.4, the basic rotating Puckon moves between minimum and maximum colatitudes. 
The Puckon is initially moving toward the North Pole, and one can see that its meridional 
velocity vanishes as it changes its direction and moves toward the South Pole. Its azimuthal 
velocity is negative throughout, i.e., it does not change its sense of rotation. Numerically, the 
Puckon moves between minimum and maximum colatitudes of (to five significant digits) 1.1033 
and 2.0384, in comparison with the corresponding values of 1.1031 and 2.0385 predicted by the 
formula given in Subsection 3.4. Also, the numerical result for the period of the Puckon's 
motion (averaged over the first five periods) agrees to six significant digits (5.60858) with the 
value given in Subsection 3.4. Thus, these numerical results for the ODE dynamics are accurate 
to between four and six significant figures. 
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Figure 4: The evolution when the initial meridional velocity is zero and the initial azimuthal 
velocity ug is a Gaussian. The larger picture shows the magnitude of the velocity. Two rotating 
Puckons have emerged at the time shown, and two more are in the process of emerging from 
the Gaussian. 
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10 20 30 10 20 30 

time time 



Figure 5: A single rotating Puckon for length scale a = 1. The initial parameters are $ = 1.5, 
P = —1, M = —2. The PDE simulation shows about half of a period of the Puckon's motion, 
and the ODE simulation shows about five periods of the Puckon's motion. The meridional and 
azimuthal velocity components and ue are shown for the PDE simulation. 
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5 Generalising to Other Surfaces 



We now wish to extend these Hamiltonian reductions for diffeons on the sphere to other surfaces. 
For the moment, we retain the rotational symmetry. 



5.1 Rotationally Symmetric Surfaces 

Any surface S with an isometric action of has rotationally invariant coordinate charts on 
which the metric is given by 

g = dr'^ + 7/;(r)2d(92 

where {r,6) E {rmin-,fmax) x (0, 27r). We may take r^m = to describe a fixed point of the 
motion. Let V be the Levi-Civita connection with respect to this metric and set 

1= (1 + A^). 

By the Bochner-Wietzenbock theorem on 1-forms we compute 

I = l + A'^-Ric 



1 + 



We also have the explicit relation, 



Suppose that is a rotationally invariant vector field on S. That is, suppose 

u = u^{r)dr + u^{r)de , 
then the singular momentum solution ansatz for EPDiff, 

lu, = ^^'' ~^h r + ^{R)5{r - R)d9 , (32) 
V'(-R) 

is a system of ordinary differential equations in r. We may solve this system to obtain two 
rotationally invariant, continuous functions G'',G^ of r and i?, within our coordinate chart, 
which are symmetric in the two variables and which satisfy 

I fGVr, R)dr + G^(r, R)de) = ^^""'J^h r + M)6(r - R)de . 

The two Green's functions C", are related by 

G\r,R) = ij{r)^{R)G'\r,R). 
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They also have a jump in the derivative along the diagonal r = R. Thus the solution of (32) 
for the velocity u is 

u = {G'' {r, R)dr + ij{r)ij{R)G'' {r, R)de)^ 

The EPDiff equations on the rotationally symmetric surface S are 

m = In, = 5^^^^--^(P,dr + M,d0) (33) 

i=l '^v-^v 

= _ + ad^m, (34) 
and these equations extremise 

^u] = ^J^{\u\l + \Vu\l)dvolg 
= {u,Iu)dvolg 
= I m{u)dvolg . 
After the Legendre transformation we arrive at the Hamiltonian 



H[m] = / m{u)dvol — £[u] 
~ 2 J "^('")'^^°^' 



We already know that the solution to (33) is 

^ G^(r,fi 

that is, 



u 



_. d_ d_ 

^ dr do 

Because m is only supported within our coordinate chart and is invariant under the action of 
S^, we have 

1 ^ 

H{m) = -2TrY, (PiMRi) + M^ueiR^)) 

4=1 
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where 27r arises as the vohime of S^. Choose a vector field w = f{r,0)-^ + g{r,9)-^. Again 
since m is supported within our chart, we know that 

^ miw)dvolg = fj (PifiRi, 0) + Mig{Ri,9)) dO. 



We also have 



d X d 

[W, U] = {fu[. - Urdrf - Uedef) ^ + {fu'g - UrdrQ 



Consequently, the EPDiff equations yield 

= / ( ^{m{w)) + m{[w,u\)\ dvol. 



-iP^fiRi,e)+MigiR^,e))d9 

2tt 



+ / (Pi (/n^ - Urdrf - uedef) + Mi {fu'g - UrdrQ - uedeg)) I^^r d6' 

^0 



2-K 





r2-K 



PifiRi, 0) + PiR^drfiR^, 6) + M^g{R^,e) + MiRidrg{Ri, 9) dO 



+ 



[ " [Pi (/< - Urdrf) + Mi [fu's - Urdrg)) d0. 

Jo 



Comparing coefficients of f,drf,gdrg now yields 

= Pi + Piu'^iRi) + M,u'g{R^ 

= PiRi - P^UriRi) 

= Mi 

= MiRi - MiUr{Ri) 

Now since 

N 

Ur{Ri) = Y,G''{RuRj) 

and 
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we see that 
d ^ 



Similarly 



This finally yields 



d 



N 



dRi 

N 

E 

j,k=i 

N 

E 

j,k=i 



j,k=i 



dRj dC^ir^R) 



dRi dr 
' dRj dG'{r,R) 



dRk dG''{r,R) 



r=R,,R=R,. 



dRi OR 



dRj dr 



+ 



dRk dG'\r,R) 



r=Rj,R=Rk 



dRj dr 



r=R^R=R. 



r=Rk,R=Rj 



dG^{r, R) 



r=Rj,R=Rk 



N 



dr 

dG''{r, R) 



+E 

r=Ri,R=Rk 



5G^(r, R) 



since G^ is symmetric in its arguments 
^ / dG^-{r,R) 

j,k=i \ 
N 

E 

k=l 
N 

^E 

2u',iRi). 



dG'^ir, R) 



dr 



r=Rk,R=Rj 



dr 



r=Ri,R=Rj 



dr 



r=R,,R=Ri 



1 d 



N 



1 d 



N 



Ri 
Mi 



2dR, 

Ur{Ri) 

0. 



(35) 

(36) 
(37) 



Thus, we have proved the following realisation of the momentum map in [5]. 



Proposition 5.1 (Canonical equations for rotationally symmetric diflfeons) The param- 
eters Pi , Ri for diffeons on a rotationally symmetric surface S with constant positive curvature 
and metric given by g = dr^+^(r)^d0^ satisfy Hamilton's canonical equations with Hamiltonian 
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given by 

F(P,R;M) = ^(^j m{u)dvol-i[u]^ 
1 ^ 

i=l 

A solution to (35,36,37) is an example of an N-diffeon on a rotationally symmetric surface S. 

We notice that the critical points of H are those R = {Ri, Rn) such that G^{Ri, Rj) = for 
ah i,j. But the value of H at these critical points is always since H depends on . 

Remark 5.2 The only 1-dimensional Lie Groups are M and , and while we have considered 
the isometric action of on a surface, we have to he more careful with translation invariance 
because the group ceases to be compact. However we can consider 6, which parameterises each 
orbit r = const, taking values in {—L,L) rather than (—00,00) for the full translation group. 
This makes the integration finite. 

5.2 Rotationally Invariant Diffeons on Hyperbolic Space 

Hyperbolic space has a richer structure than either the plane, or the sphere. Indeed, the 
isometry group of the sphere is SO (3); so all isometrics are rotations. We have already examined 
the case of rotationally invariant diffeons. (These are the Puckons.) The isometry group of 
the plane MMs k R^, the semi-direct product of rotations and translations. Translational 
invariance yields the direct product of the original 1-dimensional peakons, whereas rotational 
invariance yields the rotating circular peakons developed in [8]. However, the isometry group 
of hyperbolic space is PSL(2,M), which is not compact and contains three different types of 
isometry. These are the rotational, translational and horolational subgroups. 

The Hamiltonian for hyperbolic A^-diffeons with rotational symmetry. We have 
already done most of the work for the rotationally symmetric case. All that remains is to write 
the hyperbolic metric in a conformally flat way and from it deduce the Green's functions G^ 
and G^ . The model is familiar: we use the Poincare disc D = {{x,y) S x"^ + y"^ < 1} in 
polar coordinates with the metric 

i.e. the conformal parameter is 

^(^) = ~, 2 • 

1 — 
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The Bochner-Wietzenbock formula states that on TD 

A^^ = - 1 , 

so that 

I = 1 + A^ = 21 + A^. (38) 
We wish to find Green's functions G^' and such that 

I(G'^(r, R)dr + G''(r, R)d9) = 6{r - R) (^^dr + RdO^ . 

We also ask that these Green's functions be finite at the origin, continuous and symmetric in 
their variables. Explicitly, we have 

A'^ {a{r)dr + b{r)d9) = -^^ dr - r|- (-L|) d0. 

Upon setting 

F{r) = 2 - 2r^ + 8 log(r)r^ + 

we write 

^ y^^) = ■ / F max r,i? , 

16 mm(r, R) 

and hence 

G%r,R) = RrG''{r,R). 

Thus, we have shown the following. 

Proposition 5.3 (Hamiltonian for rotationally invariant hyperbolic A^-diffeons) 

The Hamiltonian for rotationally invariant N-diffeons on hyperbolic space is given by 

fl(P.R;M) = ^ E p(ii.y^fl.'p(Ji,y^fl, {P.P,R.R,G'(R..R,) + ^G'(R..R,)) 

N 



2 p{R,YRip{RjYRj ^ * ^ * J ^ * ^ 



The hyperbolic diffeon. Let us examine the behaviour of the basic rotationally symmetric 
hyperbolic diffeon. The Hamiltonian is 

32ii2 V ^ J 

The function F{R) /i?^ is positive and continuous on (0, 1), and it has a singularity at i? = 0. Its 
derivative is negative on (0, 1) and its image is (l,oo). Consequently, F{R)/R? is an invertible 
function (0,1) — > (l,oo). 
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Assuming M / 0, on the level set H{P, R) = we find that 

2 _ 32A'2 

So i? = 0, whenever 



F{R) 32K 



2 



i?2 M2 • 

Hence, whenever 32i^2 > m'^ > 0, precisely one turning point exists for R. This proves the 
following. 

Proposition 5.4 Rotating dijjeons with > on a rotationally symmetric hyperbolic sur- 
face do not exhibit periodic behaviour. 

For the irrotational diffeon, M = and we have 



p2 



F{R) 

Since F{R) > for all R £ [0, 1], the irrotational diffeon has no turning points of R at all. 

5.3 Horolationally Invariant Diffeons on Hyperbolic Space 

We now turn to a subtly different problem. So far we have been using solely the rotation group 
(the circle) to produce symmetric diffeons. Now we consider a different subgroup of hyperbolic 
isometrics. This time it is expedient to use the upper-half-plane model of hyperbolic geometry. 
That is, 

n = {{x,y) eR^\y>0} 

with the metric 

dx^ + dy^ 
9H = 2 ■ 

The group we consider is the horolation group {x,y) i-^ {x + a,y). This is a unique type of 
isometry in planar geometries and is, figuratively speaking, a "rotation about infinity". We 
can see immediately that the orbits of this group action are the lines y = const, and that we 
seek diffeons which are independent of x. In this situation by (38), 

I{a{y)dx + b{y)dy) = - i2ya - -^y^ - 2-Q:^y j dx + (26 - -^y^ - 2—yj dy 
Thus the solution to 

I{G{y,Y)dx) = 6{y-Y)dx 

is 

G{y,Y)- 



3 max(y, Y)"^ ' 
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and this solution also solves 

I{G{y,Y)dy) = 6{y-Y)dy. 

Next, we consider the compactness issue. Since the coordinate x ranges from — oo to oo we 
know that functions that only involve y cannot be integrated across all of TC. Thus the Diffeon 
Hamiltonian cannot exist over the whole space and the theory of section 5.1 does not apply. 
Let us restrict ourselves to the strip Bl = {{x,y)\y > 0,\x\ < L} C TC. By doing this we can 
apply similar calculations to those used in section 5.1 to Bl, but we must apply them only in 
the case of vector fields which are tangent to the boundary x = itL. 

Thus if 

N 

m = Y,Siy- Y,)Y^{Midx + Pidy) 
1=1 

then we arrive at the Hamiltonian 

N 



H{P, Y) = G{Y,,Y,)Y^Y^{P,Pj + MiMj] 



N 



LY ^'^'^Y'^' iP.P.+M.M,) 



and the condition that Mi = 0. 

For the 1-diffeon, the Hamiltonian is 



y3 

H{P,Y) = L — + m2) 
3 



so for H{P, Y) = const = K, the only turning point of Y is at 

Y = 

or at 0, if M = 0. 



3K 

Tap 



Remark 5.5 Joining the edges of Bl together yields a surface called Gabriel's Horn, whereby 
the translation invariant diffeons on Bl become rotationally invariant diffeons on Gabriel's 
horn. 



5.4 Translation Invariant Diffeons on Hyperbolic Space 

Translations on the Hyperbolic plane are characterised by having two ideal fixed points. The 
version of the metric we use in this case for the upper half plane is 

g-H = dr + cosh^ rd9^ , 
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where the subsets formed by r = const are the orbits of the translation. Again, the orbits are 
not compact, so we limit ourselves to the strip Bl = {(r, 6) £E?\r > 0, \9\ < L] dTi. 

In this situation, the Green's function G(r, R) solving 



lG{r,R)dr= ^^\lh r 
cosh(/tj 



is given by 



G(r, ^) ~ 2 tanh(min(r, R)) cosh(max(r, R)) + cosh(r) cosh(i?) arctan(e™™^'''^^). 

Due to the complicated form of the Green's function, we will not write down the explicit form 
of the Hamiltonian for the A^-diffeon. However, the Hamiltonian for the 1-diffeon is, 

HiP,R-M) = (P^cosh^R + M^) 

2 cosh R 

= ^ (tanh R sech R + 2 arctan(e^)) {P^ cosh^ R + M'^). 



— (tanh R sech R + 2 arctan(e^)) P cosh^ R , 

^ (p2 (cosh R + sinh 2R arctan(e^)) + sech^ R) . 

It can be shown that there are no equilibria for the motion unless M = 0, where the line P = 
is the critical point set for H. However, 

i7(0,fl) = 0, 

for each value of R. So the line P consists entirely of stationary points, and since away from 
this line no points exist for which i? = 0, there can be no periodic behaviour of 1-diffeons in 
this case. 



In this case 

dH 

dH 

OR 



Extending EPDifF on Hyperbolic Space Our study of EPDiff on hyperbolic spaces so far 
shows that the behaviour of hyperbolic diffeons seems less interesting than the rich dynamical 
structure of multiple Puckons interacting on the sphere. This is because hyperbolic diffeons 
will bounce only once before heading off to infinity, in the cases we have studied so far. 

However, hyperbolic space is related intimately with the study of curves of genus at least 2: 
it forms the universal cover for all such Riemann surfaces. So far, we have been unable to 
examine the case of curves of genus > 2 simply because their isometry groups are so small 
(usually discrete) that we simply cannot find a 1-parameter subgroup by which we may reduce 
the EPDiff equations to ODEs. 
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Any high genus Riemann surface can be reahsed as the quotient of the hyperboUc plane by 
a tihng, or hyperboUc lattice [17]. The complex structure of the Riemann surface depends 
upon the shape of the lattice. A surface T, of genus 7 > 2 is determined topologically by the 
symmetry group of the tiling; but there are a variety of complex structures which may be put 
upon S to turn it into a Riemann Surface. If F is the symmetry group of the tiling, then the 
space of complex structures which may be put on T, is given by the space 

F-invariant quasiconformal maps of B 

Tp = . 

Hyperbolic isometrics of B 

Here B is the Poincare disc (the unit disc in the complex plane). In the space Tp a quasicon- 
formal map is a generalisation of a conformal map in which although angles are not preserved, 
the angular dilation is uniformly bounded, see [1, 4]. 

The inclusion of a subgroup F' in F induces a natural inclusion of in Tr' , so it makes sense 
to consider the universal Teichmiiller space as the Teichmiiller space Ti associated with the 
trivial group. The universal Teichmiiller space is endowed with the Weil-Petersson metric, thus 
providing the framework essential for the theory of EPDiff. Furthermore, the restriction of the 
maps given in the definition of Teichmiiller space to the boundary of the disc provides us with 
the essential method by which we reduce the EPDiff PDEs to ODEs. This study is beyond the 
scope of this paper and we present the results in the next, [7]. 

6 EPDiff on Higher Dimensional Manifolds with Symmetry 

Our calculation for surfaces in the previous sections impose a rotational symmetry of the 
solutions. That is, the solutions are invariant under a circle action, and this effectively reduces 
I in (38) to an ordinary differential operator. As a result, we are able to find and 
which are invariant under the group action, continuous and have a jump in the derivative 
along an orbit. This motivates the study of higher dimensional Riemannian manifolds {T,,g) 
which possess an isometric action of the Lie group G such that the quotient space T,/G is a 
1-dimensional orbifold. 

6.1 Principal Bundles over 1-Dimensional Manifolds 

We consider first the case of principal bundles over 1-dimensional manifolds. It is a well known 
fact that 1-manifolds are diffeomorphic to M, [0,oo), [0,1] or S^, see for example pp55-57 of 
[13]. Intervals of M are all contractible, we know that any principal bundle over them is trivial. 
The principal i^T-bundles of are up to isomorphism in 1-1 correspondence with tti{BK), the 
fundamental group of the classifying space BK of the Lie group K. However, by the long 
exact sequence of homotopy groups (see [11]), we know that tti{BK) = 7ro(i^), thus if K is 
connected, then any Principal ET-bundle over is trivial. 

The upshot here is that if S is a Riemannian manifold with a free isometric action of K and 
T,/K is a 1-dimensional manifold A^, then T, = K x N . 
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6.2 Warped Products 



Let K he a compact semi-simple Lie group of dimension n — 1 and S — > / C M a principal 
-fC-bundle over the open interval / (in fact in what follows, we may also take / to be the circle 
S^). Suppose S has a Riemannian metric g preserved by K, and that coordinates exist such 
that the metric has the form of a warped product 

g = dr^ + ipirfgK 

where gx is a bi-invariant inner product on TK and V is a iC-invariant function on S, i.e. 
a function of r alone. The coordinate r parameterises the orbit of if in S. Let V be the 
Levi-Civita connection with respect to g. 

This situation is rich enough to produce interesting results. In what follows will be a 

gK orthonormal basis of t and the left-invariant vector field generated by ^ S We write 
Xi = X^- , these form an orthonormal frame of TK; set 6^ to be the coframe dual to the Xi . 
Let be the Levi-Civita connection of K with respect to gx- 



Proposition 6.1 We have 

A^^Xi = kXi. 



This follows from the identity on K, 

Thus, A^^^ will be the Casimir operator associated with the adjoint representation of -fC on 6 
applied to the frame {-'^^i}. The Casimir however is just a constant multiple of the identity, the 
constant k being 

^ ^ ^ _ dim3(iC) 
dim? 

where j{K) is the Lie algebra of centre of the group K. 

The geometry of warped products is reasonably well known, [2, 14]. As an exercise in Rieman- 
nian geometry, one may deduce (see pp58-59 of [2]) 



Proposition 6.2 Let f be a K -invariant function on S and suppose ip{r) = e'^^'"^. The con- 
nection Laplacian associated with g satisfies 



Qj.2 Qj. Qrf. \ Qy. 
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These two propositions show us that the X-invariant equations 

Ia°(r)dr = ^^''j^^' dr (39) 

Ia\r)ip{r)e' = 5{r - R)e' (40) 

are n ordinary differential equations in r. Thus, there are n functions G°{r, R),G'{r, R) that 
solve (39) and (40), respectively, are continuous on r = i?, and are symmetric in the sense that 

G^{r,R) = G^{R,r) 
G\r,R) = G\R,r). 

Indeed, the second identity in Proposition 6.2 shows that G* is independent of the index i. 
Hence, instead of G*,we shall write G^ . 

Now let m be the measure-valued one-form 

-^^■ = E^^J;(7j^(^^dr + M^ (41) 



The solution to l{uj)\, = is 



"J 

N 



|:(p,0«(r,7^)| + ^jJ|^G''-(,«.)X, 



Note, we are not using the summation convention here. Also denote 

for any C ^ ^- Any vector field on S will be the sum of such vector fields w. Then we may 
write the vector field commutation relation. 

Thus 

N 



l^m^{w)dvolp = f^J^l %(R^f (^'■^° + ^^^f^axiXc, X,)) ^{r)drdvolK 

N 

= V / (PifiR^, h) + Nhf{R,, h)cA dvolK 
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where C = X^IL/ C''Cj- Especially, this yields the ad* relation needed for expressing EPDiff in 
this setting, 



m^{[w,Uj]) 

= iZ f Pi f/°(^^' ^)^i'(^i) - u^i{Ri)drf{Ri, h) - Uj{Ri)X,{f){Ri, h)) dvolK 

N 

- V / Mi (u^jiRi)drfHRu h)e + u,Xjif<){R,, h)e - f{Ruh)u'ARi)) dvolK ■ 

Before we proceed, we note that while we were dealing with the circle, terms such as ugdgf 
disappeared when we integrated over the circle, as one would expect from Stokes' theorem. 
Here, however, we have the terms UjXj{f^) and UjXj{f^). One asks, "Do these vanish when 
we integrate over the group K7" As w shall see, the answer to this question is, "Yes." 

From standard Riemannian geometry, one has 

Proposition 6.3 For any compact Riemannian manifold N , and any vector field X and func- 
tion (f) on N we have 

I X{(f>)dvol= I (pdivXdvol. 

From this formula we see that, in the present situation. 
Proposition 6.4 Given any k £ t and any function (p on K. 

X^{(t))dvolK = ^ . 

K 

This relation follows because X^ generates volume preserving (metric preserving!) diffeomor- 
phisms of K. 

Hence, upon integrating over K, the terms UjXi{f^) and UjXj{f'~-) will vanish because Uj is 
JT-invariant. Thus, we have 

m^{[w,uj]) = V/ Pi(f''iRi,h)ufiR,)-u^jRi)drf'^iRuh))dvolK 

N 

- ^ (u^JiR^)^rf^iR^, h)^ " ^^(i?.)) dvolK 



i=l 
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This formula implies the EPDifF equations, 



f f dm^ \ 
= [-^[w)+m^{[w,uj])\ dvolp 

= E / ^ {Pif^R^^ M + PiRi^.fiRi^ h) + Mif<{R,, h)C^ + M,R~f^{R„ h)C^ dvolK 

N 

+ P^ (•^°(^*' - U^jiR^)^rfiR^, /l)) dvolK 

i=l ''h&K 

N 

- ^'^^ (u']{Ri)^rf^{R^, h)C - f{R^, h)u'^{Ri)) dvolK 

i=l Jh&K 

Again, we compare coefficients of f^, f^, and -^f^ to find 

= Pi + Piu']\Ri) + Miu'j{Ri), 
= PiRi - PiU^Ri) , 
= MiC' , 

= MiRiC^ - Miu]{Ri)QK 
Therefore, we have proven the following. 

Theorem 6.5 (Hamiltonian diffeon reduction for A'-invariant warped product spaces) 

The warped product diffeon parameters in equation (41) satisfy 

1 d ^ 

Pi = {PkU^Rk) + MkU,{Rk)) , (42) 

* k=l 

Ri = u%R,), (43) 
Mi = 0. (44) 

As guaranteed by the momentum map property of the diffeon solution ansatz (2) for EPDiff, 
these are Hamilton's equations with collectivized Hamiltonian Hj/{vo\K) given by 

iJj(P,R;M) = - I m^{uj)dvolp 

1 ^ 
= -vol (K) {PiU%Ri) + MiU.j{Ri)) 

4=1 

1 ^ 

= -vol {K) {Piu]{Ri) + MiU.j{Ri)) 

i=l 

Notice that Hj is independent of j; so we may write H^*" for Hj. 
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Thus, given any G t a difFeon exists that moves on S with motion described by Hamilton's 
equations with the Hamiltonian given by H. The diffeon itself is "rotating" in the direction 
determined by C, with constant angular momentum determined by M = (Mi, . . . , Mj\f). 

Remark 6.6 One may repeat the whole procedure replacing K with a compact symmetric space 
K/ H for some closed Lie subgroup H of K. The only significant changes would he that the 
9^ would become local on K/H rather than global. However, all the propositions will remain 
true because of the intimate relationship between symmetric spaces and Lie groups. Thus, for 
example, we would expect diffeon behaviour on S"" to he similar to the Puckon behaviour on 
since 

ggn = dr^ + (cos^r) c?s— 1. 

6.3 Singular fibres 

We have so far dealt with a free action of a group on a manifold such that the quotient space 
is a 1-manifold. Thinking back to the case of E = S^, we see that the action of the circle 
was not completely free, because there are precisely two points (the poles) which are fixed 
under the group. Away from these fixed points, the sphere is a principal fibration, and the 
Puckons "bounce" off the poles (provided they are not rotating). For the general manifold S, 
the situation could become much more complex. 

For example, in the situation of the previous section, we see that problems arise if we choose a 
diffeon which is rotating in the direction C € t and find that vanishes at a point p. However, 
if the diffeon is not rotating, and the vector field ^ vanishes nowhere, then all the previous 
theory holds. The theory also holds for the warped product of the line (or circle) with the flat 
metric with an Einstein space, upon using harmonic coordinate charts (p285 of [14]). 

7 Conclusions 

We provided a canonical Hamiltonian framework for exploring the solutions of the EPDiff equa- 
tions on surfaces of constant curvature. This framework used symmetry and the momentum 
map for singular solutions to reduce the EPDiff integro-partial differential equation to canon- 
ical Hamiltonian ODEs in time. We specialized to the case of the sphere and provided both 
numerical integrations and qualitative analysis of the solutions, which we called "Puckons." 

The main conclusions from our numerical study were: 

• Momentum plays a key role in the dynamics of Puckons. Radial momentum drives Puck- 
ons to collapse onto one of the poles, and angular momentum prevents this collapse from 
occuring. Puckons were found to exhibit elastic collision behavior (with its associated 
exchanges of momentum and angular momentum, but with no excitation of any internal 
degrees of freedom) just as occurs in soliton dynamics. 
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• Puckons without rotation may collapse onto one of the poles. This collapse occurs with 
bounded canonical momentum and the radial slope in velocity appears to become vertical 
at the instant of collapse. 

• For nonzero rotation, Puckon collapse onto one of the poles cannot occur and the radial 
slope in velocity never becomes infinite. 

• Head-on collisions between two Puckons may be accompanied by an apparently vertical 
radial slope in velocity which forms in finite time. 

The main theoretical questions that remain are: 

• Numerical simulations show that near vertical or vertical slope occurs at head-on collision 
between two Puckons of nearly equal height. A rigorous proof of this fact is still missing. 

• It remains to discover whether a choice of Green's function exists for which the reduced 
motion is integrable on our 2N dimensional Hamiltonian manifold of concentric rotating 
Puckons for > 1. 

• It also remains to determine the number and speeds of the rotating Puckons that emerge 
from a given initial condition. 

All of these challenging theoretical problems are beyond the scope of the present paper and we 
will leave them as potential subjects for future work. 

We applied these ideas to hyperbolic spaces, as well. This led to rather simple reduced dynamics 
with only a limited number of possible collisions. We suggested a new departure for hyperbolic 
space, based on Teichmiiller theory, which we shall investigate elsewhere. 

Finally, we answered an outstanding question by generalising the momentum map to the case 
of diffeons with {n— 1)— dimensional internal degrees of freedom by using the theory of warped 
product spaces. 

In summary, we identified and analysed cases where imposing an additional translation sym- 
metry on the solution reduced the canonical Hamiltonian dynamics of the singular solutions of 
EPDiff on Einstein surfaces from (integral) partial differential equations to Hamilton's canon- 
ical ordinary differential equations, in time. We extended our methods for surfaces to "mostly 
symmetric" manifolds in higher dimensions by using warped products. 
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